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ABSTRACT 

Previous to this report, the computer program of NASA SP-273 and NASA TM- 
86885 was capable of calculating theoretical rocket performance based only on the 
assumption of an infinite area combustion chamber (IAC). An option has been added 
to this program which now also permits the calculation of rocket performance based 
on the assumption of a finite area combustion chamber (FAC). In the FAC model, 
the combustion process in the cylindrical chamber is assumed to be adiabatic, but 
nonisentropic. This results in a stagnation pressure drop from the injector face to 
the end of the chamber and a lower calculated performance for the FAC model than 

the IAC model. 


INTRODUCTION 

The calculation of theoretical rocket performance involves a number of assump- 
tions. For the same propellant and operating conditions, theoretical performance 
can vary depending on which assumptions are used. Rocket performance calculated 
by the computer program of references 1 and 2 assumes adiabatic combustion in 
an infinite area combustion chamber (IAC) followed by isentropic expansion in the 
nozzle. In order to have a more realistic model, this supplement to references 1 
and 2 presents an additional option for calculating rocket performance based on the 
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assumption of adiabatic combustion in a finite area combustor (FAC) followed by 
isentropic expansion. Two input options are available for the FAC problem. Input 
option 1 is an assigned contraction ratio, while input option 2 is for an assigned 
mass flow per unit combustor area. 

The addition of this new FAC option required changes in only two of the sub- 
routines of the reference 1 program; namely, ROCKET and RKTOUT. A short 
description is given herein of the equations and iteration procedures. Three sample 
cases are given to facilitate the discussion of input options, output, and analysis of 
the effect on performance for the assumptions of FAC and I AC. 


FINITE AREA COMBUSTION 

Combustion in a rocket chamber is a nonisentropic, irreversible process. During 
the burning process, part of the energy released is used to raise the entropy, an 
undesirable form of energy inasmuch as it is unavailable to do work. This energy 
utilization loss is reflected in a total pressure drop as the gases are being accelerated 
from the beginning of the combustion process (at or near the injector face) to the 
end of the chamber. The combustion process may still be considered to be adiabatic; 
however, due to heat not being added at constant pressure during combustion, “the 
energy available for producing nozzle exit velocities is less than exists under ideal 
conditions of negligible chamber velocity.” (ref. 3). Calculated rocket performance 
will therefore be less for the model of FAC than for the more commonly used ideal 
model of I AC. 

A sketch of a rocket is given in figure 1. The positions in figure 1 are numbered 
in the same order as they appear in the performance output tables. The entrance 
to the finite chamber will be referred to as the injector face and will be indicated 
by ‘inj’ or ‘1’ as subscript. The end of the finite chamber (nozzle entrance) will 
be indicated with the subscript ‘4’ or ‘c\ The infinite area position is indicated by 
subscript ‘inf’ or ‘2’ while the throat is indicated by subscript ‘t’ or ‘3’. 


Equations 


Unless otherwise stated, the International System of Units (SI Units) is used. 
The relationship of forces between points 1 and 4 for the nonisentropic process of 
combustion in a finite area is given by the following equation (ref. 3, p. 81) 


P + 




( 1 ) 
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where P is pressure, A is the combustor area, rh is the mass flow rate, and u is 
velocity. Equation (1) may be written as 

(P + pu 2 )! = (P + pu 2 ) 4 (2) 

by using the continuity relationship 

m= pAu (3) 

where p is density. 

When velocity at the injector face is negligible, equations (1) and (2) reduce to 

Pi = Pin i = [ P + n y) = + P u2 h (4) 


Iteration Procedure 


An iteration procedure is required to satisfy Eq. 4. Two input options are 
available for FAC. In option 1, the contraction ratio is assigned. In option 2, 

the mass flow rate per unit combustor area ^ is assigned. The iteration procedure 
for option 1 is simpler and therefore will be described first. All of the first four 
points shown in figure 1 are involved in the iteration procedure. Thermodynamic 
parameters at point 1 are obtained by a combustion calculation (HP problem in 
reference 1). Starting with an estimated value for P 2 , calculations are then made 
for points 2, 3, and 4 (the assigned contraction ratio) as would usually be done for 
infinite area combustion, throat, and an assigned area ratio as described in reference 
1. A check is made to see if equation 4 is satisfied to within the following tolerance 



If Eq. 5 is satisfied, then the calculations for the finite area combustor are 
complete for points 1, 2, 3 and 4. Calculations are then continued if other values of 
pressure ratio and/or area ratio have been specified in the input dataset. If Eq. 5 
is not satisfied, then an improved estimate for P 2 is obtained as described in a later 
section and the procedure for points 2, 3 and 4 is repeated until Eq. 5 is satisfied. 

A similar procedure is used for option 2 as was described for option 1. However, 
the contraction ratio is not known for option 2. Therefore, the iteration procedure 
involves starting with an estimated value for as well as for P 2 and then obtain- 
ing improved estimates for both P 2 and Not surprisingly, more iterations are 
required for option 2 than for option 1 which requires improved estimates for P 2 
only. As in the case of option 1, iteration is complete when Eq. 5 is satisfied. 
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Initial Estimates 


A curve is given in figure 3-18 of reference 3 which relates -jp- with 4* for an 
assumed value of 7 = 1.2. The following empirical equation was derived by fitting 
three selected points read from the curve: 

Pi = P.n, 


1.0257 - 1.2318^ 

fit. 

1 - 1.265054* 

At 


Eq. 6 is used only to obtain an initial estimate for P 2 . 

For option 1 , the assigned value of ^ is used in Eq. 6 . For option 2 , an 
initial estimate of ^ is required (see Input Option Parameters section). This initial 
estimate is used in Eq. 6 to obtain a value for P 2 , which is then used in the following 
equation to obtain an improved initial estimate for 


_ P 2 
A t 2350 f 

Ac 


(7) 


Eq. 7 was derived by starting with the relationship for characteristic velocity c* = 
P 2 4*, multiplying both sides by A c and using an arbitrary value of c* = 2350 m/s. 
Somewhat better initial estimates for both P 2 and ^ are obtained by repeating 
several times the sequence of substituting values of ^ from Eq. 7 into Eq. 6 and 

values of P 2 from Eq. 6 into Eq. 7. If the input value of ^ is so large that Eq. 
7 calculates a value less than 1 , the program will stop the calculations and print 
out the error message “INPUT VALUE OF MDOT/A = (value) IS TOO LARGE. 
GIVES CR ESTIMATE LESS THAN 1 ”. 


Improved Estimates 

For option 1, an improved estimate for P 2 is obtained by assuming that the ratio 
of the assigned value of P, nj to the current value of P mj (obtained by means of Eq. 
4) is equal to the ratio of the final value of P 2 to the current value of P 2 . This 
assumption leads to the following equation 

( 8 ) 

The use of Eq. 8 often gives such an excellent improved estimate for P 2 that it need 
be used only once to obtain convergence (Eq. 5). 


2 ,neui 


P 

P , tnj,a 
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For option 2, an improved estimate for ^ is required in addition to the one for 
P 2 and is obtained from the following equation 

A m 

Ac _ A‘ (9) 

A t A 

a c 

Inasmuch as both P 2 and £ are changing, the iteration procedure is longer for 
option 2 than for option 1. For option 2, as well as for option 1, convergence is 
complete when Eq. 5 is satisfied. 


Input Option Parameters 

Two options are available for obtaining finite combustor area performance. In 
addition to the usual required input parameters described in ref. 1 for namelist 
&RKTINP, several additional parameters are required. For option 1 these are: 
FAC = T and AC AT = some assigned value for For option 2, the additional 

parameters are FAC = T and MA = some assigned value for Option 2 also 
requires an initial estimate of A default value of AC AT = 2 is provided in the 
program for this initial estimate! However, if desired, a different initial estimate for 
AC AT may be included in the &RKTINP data. Thus, for option 2, a value of MA 
is required in &RKTINP, whereas an estimated value for AC AT is optional. 

In FAC the PCP values are relative to the injector face pressure, whereas 

in I AC, the PCP values are relative to the infinite area chamber pressure. 
Due to this definition of PCP, the assigned values of PCP must be larger than 
p in} . Otherwise, this will give values of P e larger than Pm/, which is an impossible 

condition. For example, in table III, the value of = 1.0848. If a value of PCP 
less than this had been assigned, 1.05 for example, this would have given a value of 
P — fhi = 533170 0 _ 5077810 Pa which is more than the value of P, n / = 4914900 
Pa, an impossible condition. If impossible values of PCP are inadvertantly included 
in the input data set, these values will automatically be omitted by the program 
and the following error message printed: PRESSURE RATIO OF (value) GIVES 
PE GREATER THAN PINF. OMIT THIS POINT. 

SAMPLE PROBLEMS 

Three sample problems were selected, one for IAC (case 1) and two for FAC 
(cases 2 and 3), to illustrate some input and output features and to provide per- 
formance data for a comparison of results. The input datasets for these problems 
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are given in table I and the output is given in tables II to IV. All sample problems 
are for the same propellant, o/f and chamber pressure. The propellant is H 2 { 1) 
at 20.17A' and 0 2 (1) at 90.18A, o/f = 5.55157, and chamber pressure is 5331721 
(52.62 atm). A number of assigned pressure ratios, PCP, and supersonic area 
ratios, SUPAR (^) are common to all problems. The PCP values selected are 10, 
100, and 1000. The SUPAR values are 25, 50 and 75. In addition, the FAC case 2 
has as assigned contraction ratio = 1.58 while the FAC case 3 has an assigned 
= 1332.0. The value of j- was calculated from the results of case 2 as follows: 
from table III, in the column for ^ = 1.58, p = 2.0353 and u = J, p = 654.5. The 
product, pu = 1332, is equal to j- from the continuity relationship (Eq. 3). Case 
3, therefore, should reproduce the case 2 contraction ratio of 1.58, which indeed it 
does. Cases 2 and 3 both have FAC = T. For comparison purposes, the IAC prob- 
lem includes an assigned subsonic area ratio SUBAR = 1.58. In the FAC cases, 
the output column for the contraction ratio appears before the assigned pressure 
ratios, while in IAC, the SUBAR column appears after the assigned pressure ratio 
columns. 


Output Format 

The output format previously used for IAC has been somewhat revised to acco- 
modate FAC. These revisions are as follows: 

1. The first line of the output headings are the same for both cases and now 
read as follows: THEORETICAL ROCKET PERFORMANCE ASSUMING 
EQUILIBRIUM COMPOSITION DURING EXPANSION. The second line 
for IAC reads FROM AN INFINITE AREA COMBUSTOR, while for FAC 
it reads FROM A FINITE AREA COMBUSTOR. 

2. The line following the heading which gives chamber pressure in units of psia 
has been changed from PC to PINJ for FAC and PINF for IAC. 

3. An additional line has been added for FAC which gives either MDOT/AC = 
(value) if the input data set contains an assigned value for ^ or AC/AT = 
(value) if input contains an assigned value for contraction ratio ^ . 

4. An additional row of output has been added for FAC; namely, PINJ/P (ratio 
of pressure at the injector face to exit pressure -^r*). 

5. The next row gives the ratio of pressure at infinite chamber area to exit 
pressure. The label PC/P formerly used for IAC for this row has now been 
changed to PINF/P. 
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6. The first four columns for FAC are INJECTOR, INF CHAM, THROAT, and 
CN RATIO for conditions at the injector face, infinite area chamber, throat 
and contraction ratio. The columns for conditions at the injector face and 
contraction ratio are two additional columns which have been added for FAC 
and do not appear for I AC. 

7. When more than 13 columns of data are required, the first two columns of 
data are repeated on the second sheet of output data for IAC as before, while 
the first three columns of data are repeated for FAC. 

8. For IAC, the option remains, as before, of calculating rocket performance 
based on the assumption of either equilibrium composition, frozen composition 
or both during expansion. For FAC, only the equilibrium option is permitted. 

9. An option is provided to print intermediate output pertaining to the conver- 
gence process for or ^ . This output is obtained by setting IDBUGF 1 
in namelist &RKTINP. 


EFFECT ON PERFORMANCE 

Table II presents rocket performance data for the infinite area combustor (case 
1) while tables III and IV present similar data for the finite area combustor cases 
2 and 3. As expected, the results in tables III and IV are identical (see discussion 
in SAMPLE PROBLEMS). Table V summarizes and compares some of the data in 
tables II and III. It may be noted in table V that for the same pressure ratios, the 
area ratios and specific impulse for the case of finite area combustor are less than 
for case of infinite area combustor. This is due to a loss in total pressure during the 
non-isentropic combustion from the injector face to the end of the finite combustor. 

The term 1 - ^ represents the energy utilization loss due to this non-isentropic 

in/ . 

combustion. The energy utilization loss for this particular example (contraction 
ratio equal 1.58) is about 3.12% at a pressure ratio of 10 and about 0.62% at a 
pressure ratio of 1000. There are two general trends in energy utilization losses. 
The first trend, which was just illustrated, is that for the same contraction ratio, 
energy utilization losses decrease with increasing pressure ratios p } . The second 
trend, for which data are not given in this report, is that for the same pressure ratio 
energy losses decrease with increasing contraction ratios. 

The previous numerical comparisons of table V data are for the same pressuie 
ratios for IAC and FAC. However, the area ratios are not the same. When the 
comparison in energy utilization loss is for the assumption of the same area ratios, 
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the losses are negligible. For example, as may be seen in table V, for area ratios of 
25, 50 and 75, the energy utilization loss is only 0.05% or less. 


CONCLUDING REMARKS 

Previous to this report, the computer program of ref. 1 permitted calculation of 
theoretical rocket performance based on combustion in an infinite area combustor. 
An option has now been added to this program that permits performance calcula- 
tions based on the assumption of a finite area combustor. Calculations were made 
for a typical example (If 2 — 0 2 propellant, contraction ratio of 1.58) based on the 
two assumptions of finite and infinite area combustion chambers in order to assess 
the size of energy utilization losses due to the nonisentropic combustion process 
in the finite area combustor. The comparison of an energy utilization loss term 
involving specific impulse was made at several assigned pressure ratios and several 
assigned area ratios. The comparison showed energy utilization losses of 0.6% to 
3.0% for assigned pressure ratios of 1000 to 10 respectively, whereas for assigned 
area ratios of 25 to 75, the energy utilization losses were trivial (0.05% or less). 

Further information on the code can be obtained from the authors. Contact 
COSMIC, The University of Georgia, Athens, Ga. 30602, concerning the availability 
of this program. 
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APPENDIX - SYMBOLS 


A area, m 2 


ACAT 

symbol in program for con- 
traction ratio 

Ai 

cont raction ratio (ratio of fi- 

t 

nite chamber area to throat 
area) 

At 

A. 

ratio of nozzle exit area to 


throat area 

C * 

characteristic velocity, ^ 

FAC 

finite area combustor 

HP 

assigned enthalpy and pres- 
sure problem (combustion 
at constant pressure) 

IAC 

infinite area combustor 

h P 

specific impulse with 

exit and ambient pressures 
equal. , ^ or — 

‘ kg j sec sec 


MA 


m 


symbol in program for ratio 
of mass flow rate to cham- 
ber area, 


Ac 




mass flow rate. - £ 


pressure, 


N 


PCP symbol in program for ratio 

of chamber pressure to exit 
pressure (For FAC, PCP = 
For IAC, PCP = 

SUPAR symbol in program for su- 
personic area ratio 
u velocity, ^ 

7 ratio of specific heats 

p density, ^ 


a assigned 
c combustor 
e exit 
f finite 

i infinite or ideal 


Subscripts 

inf infinite 
inj injector 

o/f oxidant-to-fuel mass ratio 
t throat 
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TABLE I. - INPUT FOR SAMPLE CASES 


Case 1 Input . - Infinite Area Combustor 


REACTANTS 
H 2. 

0 2 . 

NAMELISTS ^ m 

CINPT2 KASE= 1 , RKT = T , P=52 . 62 , OF- T , MIX- 5 .55157, 

SIUNIT=T LEND 

CRKTINP SUBAR=1.58, 

PCP= 1 0 ,100,1000, SUPAR=25 ,50,75 LEND 


100 . 
100 . 


-2154 . 
-3102 . 


L 

L 


20.17 F 

90.18 0 


Case 2 Input. - Finite Area Combustor, Option 1 

100 . -2154 . L 20 . 17 F 

100 . -3102 . L 90.180 


REACTANTS 
H 2. 

0 2 . 


NAMELISTS _ ^ m , ,, K7 

EIMPT2 KASE=2 , RKT=T , P=52 . 6 2 , OF=T ,MIX = 5 .55157 , 

SIUNIT=T LEND 

CRKTINP FAC=T , ACAT= 1.58, 
PCP=10,100,1000,SUPAR=25,50,75 LEND 


Case 3 Input. - Finite Area Combustor, Option 2 


REACTANTS 
H 2. 

0 2 . 

NAMELISTS _ m 

CINPT2 KASE=3 , RKT=T ,P=52.62,OF-T, MIX- 5 . 55157 , 

SIUNIT=T LEND 

CRKTINP FAC=T , MA= 1 3 32 . 0 , 

PCP= 1 0,100,1000, SUPAR = 2 5 ,50,75 LEND 


100 . 
100 . 


-2154 . 
-3102. 


L 

L 


20.17 F 

90.18 O 
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TABLE II. - THEORETICAL ROCKET PERFORMANCE ASSUMING EQUILIBRIUM COMPOSITION DURING EXPANSION 
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TABLE III. - THEORETICAL ROCKET PERFORMANCE ASSUMING EQUILIBRIUM COMPOSITION DURING EXPANSION 

FROM FINITE AREA COMBUSTOR 
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TABLE IV. - THEORETICAL ROCKET PERFORMANCE ASSUMING EQUILIBRIUM COMPOSITION DURING EXPANSION 
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TABLE V. COMPARISON OF SPECIFIC IMPULSE OBTAINED 
UNDER ASSUMPTION OF EXPANSION FROM FINITE AND 
INFINITE COMBUSTION CHAMBERS 


Infinite 

chamber area 

Finite 

chamber 

area 

Energy Utilization loss 
due to finite chamber 

P*nf 

J 

A' 

Pin] 

r 

A± 

1 - 5 

Pe 

l >p 

(— ) 

At 

Pe 

{—) 
v sec ' 

At 


10.00 

2972.5 

2.3469 

10.00 

2925.8 

2.2253 

0.0312 

100.00 

3886.3 

12.179 

100.00 

3861.9 

11.482 

0.0125 

1000.00 

4383.2 

68.360 

1000.00 

4369.5 

64.394 

0.0062 

262.17 

4129.4 

25. 

283.25 

4128.3 

25. 

0.0005 

659.57 

4313.1 

50. 

714.18 

4312.2 

50. 

0.0004 

1131.38 

4402.6 

75. 

1225.05 

4401.8 

75. 

0.0004 


15 




INFINITE AREA CHAMBER 

© 



FIGURE 1. - THRUST CHAMBER SCHEMATIC WITH POSITIONS LABELED AS THEY APPEAR IN PROGRAM 
OUTPUT. 
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